A stochastic model of cascades in 2D turbulence 



Peter D. Ditlevsen and Jes Ravnb0l 

Centre for Ice and Climate, The Niels Bohr Institute, University of Copenhagen 

(Dated: January 4, 2012) 



O 

(N 
O 

Q 

o 

m 

Q 

U 
d 



> 

O 



o 



The dual cascade of energy and enstrophy in 2D 
turbulence cannot easily be understood in terms 
of an analog to the Richardson-Kolmogorov sce- 
nario describing the energy cascade in 3D tur- 
bulence. The coherent up- and downscale fluxes 
points to non-locality of interactions in spectral 
space, and thus the specific spatial structure of 
the flow could be important. Shell models, which 
lack spacial structure and have only local interac- 
tions in spectral space, indeed fail in reproducing 
the correct scaling for the inverse cascade of en- 
ergy. In order to exclude the possibility that non- 
locality of interactions in spectral space is crucial 
for the dual cascade, we introduce a stochastic 
spectral model of the cascades which is local in 
spectral space and which shows the correct scal- 
ing for both the direct enstrophy - and the inverse 
energy cascade. 

The scaling relations for the energy spectrum in turbu- 
lence are consequences of the energy cascade in 3D turbu- 
lence [J and the dual cascades of energy and enstrophy 
in 2D turbulence @. The cascades follow from inviscid 
conservation of energy and enstrophy (in 2D), separation 
between forcing and dissipation scales and the sweeping 
of smaller eddies by larger scale flow, showing that inter- 
actions are local in spectral space. The last assumption 
is based on the Richardson picture of energy flowing in 
a self-similar manner from large to small scales Q, and 
thus depends somewhat on the actual physical structure 
of the flow. In fact, experiments Q and simulations @ 
show that fully developed 3D turbulence has a filamented 
structure rather than a structure with eddies, or vortices, 
with smaller eddies inside as envisaged by Richardson. 

In 2D turbulence the situation is quite different : Since 
enstropy is also an inviscid invariant, the energy cascade 
is from smaller to larger scales. This inverse energy cas- 
cade is seen in both experiments @ and simulations 
The scenario fits poorly into the Richardson picture. The 
2D turbulent flow is, besides the lower dimensionality 
also different in its spatial structures from 3D turbulent 
flow. A 2D turbulence analog for the Richardson pic- 
ture could be as follows: The flow is characterized by 
well localized energy-containing vortices where the flow 
in between the vortices is characterized by a strong shear 
accounting for the enstrophy dissipation. Vortices can 
merge, leading to even larger scale structures. For this 
picture the dual cascade can be explained in a simple 
heuristic scenario: Consider a vortex of linear scale R ro- 
tating as a rigid body with rotational speed CI. Thus the 
velocity is Uj(r) = eijiCljri for r < R, falling off rapidly 
for r > R. Cli is perpendicular to the plane of the flow. 



The energy of the vortex is E = (ir / 4)S7 2 R 4 and the en- 
trophy is Z = 4Trn 2 R 2 . 

Consider now a flow of two such vortices of linear size 
R. Assume that they scatter in process after which two 
vortices of linear sizes, say, R/2 and 2R emerge. This is 
schematically shown in figure [TJ These new vortices have 




FIG. 1: Schematic cartoon of the scattering of two vortices of 
radii R in a 2D flow. The lower vortex is stretched in the flow 
of the upper vortex. This results in growth to size 2R of the 
upper vortex while a smaller vortex of size R/2 is scattered of. 
If the vortices are considered to perform rigid body rotations, 
the big upper vortex contains most of the energy while the 
small lower vortex contains most of the enstrophy. In this 
way the energy is cascaded to large scales while enstrophy is 
cascaded to small scales. 

rotational speeds fii and O2 respectively. From energy 
and enstrophy conservation Ox and O2 are determined, 



£7(7r/4) = 2V 2 R 4 = nj(R/2) 4 + n 2 2 {2R) 4 
Z/(4tt) = 2fl 2 R 2 = Ql(R/2) 2 + nl(2R) 2 
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from which we get ilj = 32fl 2 /5 and Sl\ = O 2 /10. The 
energy is then redistributed such that E\ = E/5 and 
E2 = 4_E/5, while the enstrophy is distributed such that 
Z\ = 4Z/5 and Z 2 = Z/5. Thus the energy has moved 
to larger scales while the enstrophy has moved to smaller 
scales. The role of coherent structures and scattering of 
vortices in the cascade process is not at present clear and 
this is just one of several conceptual pictures proposed 

A complementary approach to modeling the cascade 
process is taken in the one dimensional shell models intro- 
duced by Obukhov [ll| and Gledzer [l2[. In these mod- 
els there are no meaningful representation of the spatial 
structure of the flow. The flow is represented by a set of 
generalized spectral velocity components u n , associated 
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with a wave number k n = A", where A is a spectral shell 
spacing (typically A = 2). The velocity u n can be inter- 
preted as some average representation of all spectral fluid 
velocity components u(k) within a shell < \k\ < k n , 
thus the name 'shell model'. The dynamics of the shell 
models [l3|, E3] are, except for the tensorial structure, 
similar to the spectral Navier-Stokes equation: 



Un = ik n \ U* n+1 U n+ 2 - -1t*_iUn+l j7j-Un-2Un-l j 

- {vk 2 n + vik~ 2 )u n + f8 n ,n (2) 

The terms in the first parenthesis correspond to the non- 
linear advection and pressure gradient terms. The next 
term is the viscous dissipation and in the 2D case the 
large scale drag. The last term is a forcing term localized 
at some wave number k no . The shell models have two 
quadratic inviscid invariants, energy: 

n 

and a second invariant determined by the free parameter 
e: 

£ 2 =£(e-l)-> n | 2 =£fc> n | 2 , (4) 

n n 

where the last equality defines the exponent a = 
-log(e- l)/logA. For < e < 1 the factor (e - 1)-" = 
(— l)"|e — l|~ n has alternating signs for even and odd 
shell number n corresponding to a generalized helicity 
[Til Il6l | . For these parameter values the shell models are 
denoted 3D-like. For 1 < e < 2 the second invariant is 
always positive corresponding to a generalized enstrophy, 
and the models are denoted 2D-like. The enstrophy has 
the same dimension as in the real 2D flow for a = 2 and 
thus e = 5/4. The velocity in the shell models have no 
meaningful spatial structure, but the 3D-like shell models 
do exhibit a forward energy cascade, with a Kolmogorov 
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scaling (K41) relation (\u n \) ~ k n . Recent interest 
in the 3D-like shell models has been on the numerical 
finding that not only do the models show K41 scaling 
relations, the models also show intcrmittcncy corrections 
to K41 leading to anomalous scaling relations similar to 
what is seen in high Reynolds number 3D turbulence 
[T7I [l8j . From inviscid energy conservation, there is an 
exact scaling relation for the non-linear flux of energy 
(n„) = fc n A„ + i - fc„_i(e - 1)A„ = e, where e is the 
mean energy dissipation, and A„ = (Im(u n -iu n u n +i)) 
is a specific third order structure function. This corre- 
sponds to the 4/5th law of homogeneous and isotropic 
3D turbulence. 

For the rest of this paper we shall focus on the 2D-like 
models and denote E2 = Z . This case is more tricky: 
From classical scaling arguments we get the (constant) 
mean non-linear flux of enstrophy through the inertial 



range as (Ilf ) ~ *# +1 <|tt»|) 3 => (\u n \) ~ fc" (Q+1)/3 . This 
is the corresponding Kolmogorov-Kraichnan scaling for 
the shell models. Obviously, in this case, as in 2D turbu- 
lence, the inviscid enstrophy conservation also lead to an 
exact sca ling relation for a specific third order structure 
function |l0| . One heuristic argument for the transfer 
of enstrophy to smaller scales (larger wave numbers) is 
that the wave-wave interactions will tend to distribute 
enstrophy evenly over the degrees of freedom of the sys- 
tem, which is the maximum entropy state. This state 
of equipartition of enstrophy defines a different scaling 

relation; k"(\u n \ 2 ) ~ const. =>■ (|n n |) ~ k n a ^ 2 ■ 

Now, for the dimensionally correct enstrophy (a = 2) 
the two scalings are the same, so a cascade and a diffu- 
sive transport of enstrophy in quasi equilibrium cannot 
be distinguished . This is an artifact of the shell mod- 
els not present in 2D turbulence, where the spectral slope 
for enstrophy cascade is (|u(fc)|) ~ fc _1 , while for equipar- 
tition it is (|u(fc)|) ~ fc -1 / 2 . If the exponent a is differ- 
ent from 2, the scalings corresponding to enstrophy cas- 
cade and equipartition are different. For < a < 2 the 
2D-like shell models show a forward enstrophy cascade, 
while for a > 2 they show an equipartitioning of enstro- 
phy. This numerical finding could be related to how the 
typical eddy turnover time r„ depends on wave number. 
The typical eddy turnover time is simply defin ed by di- 
mensional counting: t„ = (fc„(|u„|)) _1 = (k n ^ (E n )) -1 . 
Assuming the scaling E n ~ k^, the scaling for the eddy 
turnover time becomes r„ ~ = fc r / 7+2 ^ 2 . The scal- 
ing exponents in the energy range and for three values 
of a in enstrophy range in the two cases of cascade or 
equipartition and are summarized in Table 1. 
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TABLE I: Scaling exponents 7 and k for the energy spectrum, 
E n ~ kZ, and the eddy turnover time, r n ~ k n in the cases of 
cascade or equipartition of energy or enstrophy respectively. 

In the case of a cascade (0 < a < 2) the typical 
turnover time decrease with increasing wave number, 
while in the case of equipartition (a > 2) it increase 
with wave number, leaving time for upscale (from large 
to small wave numbers) transport of enstrophy to equi- 
librate 0, [2l| . For the same reason the 2D-like shell 
models fail in simulating the inverse cascade of energy: 
For a spectrum corresponding to energy cascade the eddy 
turnover time decrease with increasing wave number (see 
Table 1), which makes the transport of energy diffusive, 
preventing the classical inverse energy cascade (indepen- 
dent of a). The situation is summarized in figure[5J where 
the energy spectra for the three cases, a = 1,2,3 are 
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shown: For a = 1 the model has a cascade spectrum 
(dashed line), for a = 3 it has an equilibrium spectrum 
(full line), while for a = 2 the two spectra coincide. 
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FIG. 2: Shell model energy spectra for a — 1,2,3, (e = 
1 + \~ a ). The pumping scale is at wave number no = 15. 
Other parameters of the simulations are: v = 10 _14 ,^i = 
100, / = 1 + i. The dashed lines are the scaling relations 
corresponding to (forward) cascade of enstrophy , while the 
full lines corresponds to equipartition of enstrophy. When the 
slope corresponding to equipartition is steeper than the slope 
corresponding to cascade (a = 3) , the model shows equiparti- 
tion, while in the opposite case (a = 1) it shows cascade. By 
the same token there is equipartition of energy in the inverse 
cascade range. 



It thus seems that one dimensional models are unable 
to generate the dual cascade phenomenon characteristic 
for 2D turbulence, which suggests that the specific spa- 
tial structure of the flow is essential for the inverse energy 
cascade. In order to investigate this or if the dual cas- 
cade is related to specific scale dependence of turnover 
times we construct a stochastic Markov chain model of 
the cascade process. As for the shell models wc define 
a chain of exponentially growing scales in wave number 
space k„ = A". The dynamical variable associated with 
each scale is the energy E n . The enstrophy Z„ is related 
to the energy as Z n = k"E n . The stochastic dynamical 
equation for E n is: 

dE n = { q n+ iE n+ i + (e — l)q n _iE n -i — eq n E n 
+ q n +iE n+ i + (e — l)g n _i.E n _i — eq n E n 
- {vk 2 n + v x k~ 2 )E n + fS n . no }dt (5) 

where E n = v / E n ^ 1 E n+ i and {(q„,q n ), n = 1,...,N} is 
a set of 2N stochastic variables: 



1/St (— 1/ St) with probability P n 
with probability (1 



where P n = min(l, <5t/t„), t„ = l/{k n y/E^) is defined as 
a dynamical eddy turnover time and St is a time interval 
smaller than the smallest time scale in the system. It is 
straight forward to verify that energy and enstrophy are 
conserved in the unforced and inviscid case. The case 
q n = 1 /St corresponds to a triad interaction where en- 
ergy is transferred from shell n to shells n — 1 and n + 1. 
The case q n = —I/St corresponds to a triade interaction 
where energy is transferred from shells n — 1 and n + 1 to 
shell n. The choice of E n as the geometric mean of the 
energies of the neighboring shells ensures that energies 
remain positive. Furthermore, in the case that the en- 
ergy follows a perfect scaling relation, E n = E (k n /ko) 1 , 
we have E n = E n and the model has detailed balance 
in the sense that a (positive) energy /enstrophy transfer 
from shell n to the two neighboring shells has the same 
probability MS M transfer in the opposite direction. 

The stochastic model energy spectra for the three cases 
a = 1 , 2, 3 are shown in figure [31 In all three cases the 
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FIG. 3: Stochastic model energy spectra for a = 1, 2, 3, (e = 
1 + A~ Q ). The pumping scale is at wave number no = 15. 
Other parameters of the simulations are: v — 10~ 17 ,!A = 
0.5,/ = 0.1. The lines for n > no are the scaling relations 
corresponding to (forward) cascade of enstrophy, while the 
lines for n < no correspond to the K41 scaling for the inverse 
cascade of energy. 



two scaling regimes of inverse cascade of energy and for- 
ward cascade of enstrophy are observed. The stochastic 
model thus, in contrast to the shell models, shows the 
same behavior of dual cascade as in 2D turbulence, so 
even though the eddy turnover time in the spectral range 
of inverse cascade of energy decrease with wave number, 
the system will not equilibrate. In the spectral range of 
inverse energy cascade where large scale energy dissipa- 
tion and smaller scale energy pumping is well separated, 
there is a statistical steady state, (n„) = e, where again 
e is the mean energy dissipation. Correspondingly in the 
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range of forward enstrophy cascade there is a statistical 
steady state, (ILf ) = rj, where rj is the mean enstrophy 
dissipation. Similar to the shell model, the mean non- 
linear transfer of energy (II n ) and enstrophy from 
shells m < n to shells m > n are easily calculated: 

(n n ) = (q n+1 E n+1 ) + (1 - e)(q n E n ) 

+(q n+1 E n+1 } + (l-e)(q n E n ) (7) 

and 

(nf) = k%({q n+1 E n+1 }-{q n E n ) + (q n+1 E n+1 }-(q n E n }). 

(8) 

Each of the terms on the right hand sides has the form 
(q n E„) = (P n E n ) = k n (En /2 ) = k n A„ or (q n E n ) = 
-{Pjtn) = -k n ({E n - 1 E n E n+1 ) 1 /' 2 ) = -k n A n and 
equations and |S} can be rewritten 

(n n ) = fc„{A(A„ +1 - A„ +1 ) + (1 - e)(A„ - A n )} (9) 

and 

(nf > = K +1 {X(A n+1 - A„ +1 ) - (A„ - A„)}. (10) 

An exact scaling relation A„ = A„ = ckfC 1 ^ 2 would imply 
(II n ) = (11^) = violating the non-zero inverse energy 
and forward enstrophy cascades. Numerical inspection 
shows that (A„ — A„)/A„ w 0.02 independent of n (and 
a). Thus we may assume a K41 scaling relation (A„ — 
A„) = Ckf7 /2 . In the range of inverse energy cascade, 
k n < k n} (or forward enstrophy cascade, k n > k nf ), 
(II n ) = e (or 0) and (nf ) = (or rj), equations (|9]) and 
(fT0]l implies: 

(n„) = Ckl+ 3 ^ 2 {X 1+3 ^ 2 + (1 - e)} = e (or 0) (11) 
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curly brackets. 
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